Collision of one-dimensional fermion clusters 
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We study cluster-cluster collisions in one-dimensional Fermi systems with particular emphasis on 
the non-trivial quantum effects of the collision dynamics. We adopt the Fermi-Hubbard model and 
the time-dependent density matrix renormalization group method to simulate collision dynamics 
between two fermion clusters of different spin states with contact interaction. It is elucidated 
that the quantum effects become extremely strong with the interaction strength, leading to the 
transmittance much more enhanced than expected from semiclassical approximation. We propose 
a concise model based on one-to-one collisions, which unveils the origin of the quantum effects and 
also explains the overall properties of the simulation results clearly. Our concise model can quite 
widely describe the one-dimensional collision dynamics with contact interaction. Some potential 
applications, such as repeated collisions, are addressed. 



Recently non-equilibrium dynamics of cold atoms has 
attracted much attention, because cold atom systems 
are ideal as isolated quantum systems which can be ex- 
perimentally designed [l]. In cold atom systems, Fesh- 
bach resonance [2| has enabled experimentalists to mod- 
ify strength and sign of interaction between atoms. Also, 
the dynamics of quantum quench ^ has been explored 
by suddenly changing the trap potential and the interac- 
tion. 

Of these experiments, the dynamics of collision and 
mixing of two fermion clouds released from spin- 
dependent traps 'i'l motivates us to study collision dy- 
namics between two fermion clusters in one-dimensional 
Fermi systems. For such one-dimensional collision dy- 
namics, a well-defined semiclassical model exists [El], so 
the fully quantum results and the semiclassical results 
can be theoretically compared. Thus the differences be- 
tween the two results, which we call the quantum ef- 
fects, can be calculated. The one-dimensional collision 
dynamics 043 has been simulated from the interest in 
that experiment. It has been shown that the sign of the 
interaction does not affect the dynamics, and that the 
non-trivial quantum effects are dominant when the in- 
teraction is strong. 

However, the origin of the quantum effects has not 
been clarified within our understanding. The complete 
analysis in the quantum effects of the collision dynam- 
ics is essential to an understanding of the quantum non- 
equilibrium dynamics, since the collision dynamics is one 
of the basic concepts of dynamics. 

In this study we simulate the cluster-cluster colli- 
sion dynamics in one-dimensional spin- 1/2 Fermi sys- 
tems with contact interaction, and calculate the quan- 
tum many-body effects in the collision dynamics. Then 
we propose a quantum collision model which contains 
only the one-to-one collision parameters. This model, 
the independent collision model (ICM), reproduces well 
the simulation results, and fully elucidates the non-trivial 
quantum effects of the collision dynamics, which ex- 
tremely enhance the transmittance at strong interaction. 



The ICM can quite widely explain the one-dimensional 
collision dynamics with contact interaction. We address 
the two applications of the ICM, interaction sign effects 
and repeated collisions. 

We conduct simulations of collision dynamics, and 
calculate reflectance and transmittance of the clusters 
in one-dimensional spin-1/2 Fermi systems. Initially n 
fermions per spin are trapped by spin-dependent poten- 
tials separately (FigH]) (n < 6 for numerically exact re- 
sults). Both trap potentials are harmonic, and they have 
the same shape but are spatially separated. The mass of a 
fermion is m, and the trap frequency of the harmonic po- 
tentials is oj (the oscillation cycle is T = 2tt/uj), and the 
interaction strength between the fermions is zero. So the 
typical width of the particle density tails is 77 = ^Jh/muj. 
We set the distance between the centers of the initial po- 
tentials as 2D = IO77 ^ 77, so that the overlap between 
the two clusters is negligible. The center of the spin-down 
(up) trap potential is x — ~D — —^j] {x — D — Sry). 
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FIG. 1. (Color online) Initial particle density and trap poten- 
tial at n = 6, and initial spin-down density at n = 1, n = 2, 
n = 3. 



At i = we suddenly change the trap potentials into a 
new shared potential, V{x) ~ imw^x^. Simultaneously 
we switch on the contact interaction between spin-up and 
spin-down fermions as u6{xd — Xu), where Xd (xu) is loca- 
tion of a spin-down (up) particle. Then the two clusters 
start moving towards each other without significantly 
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changing their shapes, and they coUide around a; = 
at t = T/4 with the average momentum |p| = muD. Fi- 
nally, depending on the interaction strength u, they are 
reflected to the initial location, or travel to the oppo- 
site location hy t — T/2, and the particle density around 
X = is almost zero again. So the number of reflected 
particles N^^u) and the number of transmitting parti- 
cles N^^{u) are obtained by counting the spin-down (up) 
particle number in x < and x > {x > and a; < 0). 
We obtain the reflectance i?n(u) — N^^{u)/n and the 
transmittance T„(m) = Nl^''{u)/n {Rn{u) -f r„(u) = 1). 
Clearly i?„(0) = 0,T„(0) = 1 and i?„(cx)) = l,r„(oo) = 
because the system is one-dimensional. 

We discretize the system to adopt the one-dimensional 
Fermi-Hubbard model, and apply the time-dependent 
density matrix renormalization group (t-DMRG) method 
to simulate the dynamics. We take 199 sites num- 
bered —99, —98, • • • , +98, +99 at regular intervals; the 
site —50 (+50) is the initial location of the potential cen- 
ter for spin-down (up) atoms. The lattice constant is 
5x — 2i5/100 = O.I77, which is small enough to neglect 
the umklapp scattering. The trap potential of the site i 
is Vi,cr{t) — ^mijp{xi — aDY (c = — for spin-down, and 
(7 = + for spin- up). The discretized Hamiltonian is 



2mSx^ 



We calculate the time evolution of this Hamiltonian start- 
ing from the ground state of the system by t-DMRG up to 
t = T/2, where the time step is 10~^T and the maximum 
eigenvalue of the density matrix discarded is £ < 10~^^. 
The simulation is conducted in the following range of 
parameters: the fermion number n < 6, the contact in- 
teraction strength < u/uc < 2^ {uc — 2fvp/m). 

Figure [2Ka) shows the reflectance Rn{u) obtained by 
the DMRG simulation for 2"^ < u/uc < 1. The fig- 



ure implies Rn{u) 



in the small u limit. We also 



plot the ratio Rn{u)/Ri{u) in Figl2jb) to evaluate the 
many-body effects. It is observed that i?„(u) ~ nRi{u) 
is approached in the limit of w — > 0. In this limit, al- 
most all particles transmit, and a particle collides with 
the n particles of different spin until it reaches the oppo- 
site side. So in the semiclassical picture, the reflectance 
is nRi{u) because the number of the reflected particles 
is approximately n'^Ri{u). Therefore there is no quan- 
tum effect in this limit 5]. On the other hand, FigHl^a) 
shows the transmittance T„(m) from the DMRG simula- 
tion for 1 < u/uc < 2^. The figure shows Tn{u) oc u^"^ 
in the large u limit. We plot the ratio Tn{u)/Ti{u) in 
FigMb), which illustrates Tn{u) ~ nTi(u) as it — )■ 00. In 
this limit, since almost all particles are reflected, in most 
cases a particle collides with another particle of differ- 
ent spin just once during a cluster-cluster collision. So 



in the semiclassical case, the number of the transmitting 
particles is almost nTi{u), so the transmittance is Ti(u). 
Therefore there are strong quantum effects in this limit 
, and the transmittance in the quantum case is n times 
larger than in the semiclassical case. 
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FIG. 2. (Color online) (a) Reflectance obtained by DMRG 
simulation Rn{u) (dots) and that by ICM R^^^{u) (lines) 
(log scale on x and y axes), (b) Reflectance ratio obtained 
by DMRG simulation R„{u)/Ri{u) (dots) and that by ICM 
R^n^{u)/RT^{u) (lines) (log scale on x axis). 



We propose the independent collision model (ICM) to 
explain the simulation results. We start from the sim- 
plest case, the one-to-one collision dynamics, which can 
be easily calculated. The initial wavef unction is |i t)i 
which means that a spin-down (up) particle is in the left 
(right). When the spin-down particle with momentum 
p and the spin-up particle with momentum — p collide, 
lit) splits into p I i t ) (reflection term) +r 1 1 i ) (trans- 
mission term) and the momenta are reversed, in which 
p = u/(iuc — u), T = iuc/{iuc — u). The theoretical 
reflectance R\'~^^{u) and transmittance Tj'^^''(u) are cal- 
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FIG. 3. (Color online) (a) Transmittance obtained by DMRG 
simulation T„{u) (dots) and that by ICM r^°^(w) (lines) (log 
scale on x and y axes), (b) Transmittance ratio obtained 
by DMRG simulation Tn{u)/Ti{u) (dots) and that by ICM 



Tl^^{u)/Tt'-''\u) (lines) (log scale on x axis). 



culated by the coefRcients of || t ) and 1 1 i ) as 



Next we regard the multi-particle cases as a se- 
ries of one-to-one collisions. In this system, the ini- 
tial wavefunction can be expressed as a Hartree-Fock 
form. We set up the single-particle wavefunctions of 
spin-down particles as ip-i{x), (p-2{x), ■ ■ ■ ,(p-n{x), and 
those of spin- up particles as (pi{x), ip2{x), ■ ■ ■ ,(pnix). So 
the initial wavefunction of the whole system is ip = 
(nl)~^ \ip^i{x^j)\ ^ \ipi{xj)\^, where the variables in the 
determinant | • | ( | ' | -|- ) belong to spin-down (up) . We 
can take these single-particle wavefunctions not only as 
the eigenstates of the initial Hamiltonian, but also as the 
localized wavefunctions created by a unitary basis trans- 
formation. It is convenient to use the localized single- 



particle wavefunctions to analyze the cluster-cluster col- 
lision. The important property of the contact interac- 
tion is that there is no difference of the shape of the 
single-particle wavefunction before and after a one-to- 
one collision but the whole wavefunction splits into the 
unchanged part (transmission term) and a spin-flipped 
part (reflection term). Therefore, if the one-to-one 
collisions occur independently, the time evolution of the 
system is described by (i) the time evolution of the single- 
particle wavefunctions and (ii) the splits of wavefunctions 
at the one-to-one collisions. The time evolution (i) oc- 
curs in each ip±i{x) and the form of ip±i{x) is changed, 
but this evolution does not change the expression of 
ip — {nl)~^ \(p_i{x-j)\ ^ \Lpi{xj)\^. On the other hand, 
at the time evolution (ii), the wavefunction splits into 
the transmission term and the reflection term. When 
the particle a and the particle b collide, ip splits into 
Tip + pF^^ (Fj" is the flip operator of and ipb)- Even 
when a and h have the same spin, that relation holds 
because r -I- pF^^ = 1 for the particles of the same spin. 

We assume that ipi{x) is localized at Xi (a;_„ < • • • < 
X-i and xi < ■ ■ ■ < Xn) and the time evolution (i) 
is expressed as the motion of Xi, and that tpi indepen- 
dently collides only with the wavefunctions of the re- 
versed momentum at the same location. This assump- 
tion, the independent ordered collision, is supported nu- 
merically later. Since the movement of Xi is similar to 
the case of the classical dynamics, the ordering of the 
collisions is the same (Fig. |4l which is discussed later). 
The initial wavefunction of the whole system is expressed 
as -0 = Ii"" • ■ ■ i~^t^ ■ ■ • t"), which means that there 
are spin-down particles at Xi {i < 0) and spin-up at xj 
{j > 0) and the ordering of wavefunction in Slater deter- 
minant is adjusted to the one-dimensional ordering of x^. 
The first one-to-one collision occurs between x-i and xi. 
So the wavefunctions splits into 

ri-.-rVHY •••> + pi---r'rUY •••> , 

but the order of x_i and xi is reversed, so the expression 



p\---r''i'r'f ■■■) + r\- 



corresponds to the spatial spin distribution. Simultane- 
ous collisions are commutative, so we can assume that 
one of them occurs earlier. The next collision is assumed 
to be between x-2 and xi (simultaneously X-i and X2), 
■) is not changed, but for the sign 
the spatial spin distribution ex- 



then \--- i-^i^t'^f ■ 
of Slater determinant 
pression is 

|...|-2|l^-1^2 . 



By omitting the ordering of Xi from this expression, we 
obtain the simplified expression of the system wavefunc- 
tion. It is possible to calculate the outcome of the cluster- 
cluster collision by calculating the system wavefunction 
after one-to-one collisions between all possible com- 
binations of Xi (i < 0) and Xj (j > 0), because the 
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single-particle wavefunctions evolve as that of free parti- 
cles, and they collide n times during the half cycle of the 
oscillation. 




FIG. 4. (Color online) Collision order in ICM at n = 3. S-ij 
is the one-to-one collision between ip-i and tpj, and tpk is the 
wavefunctions after the kth set of collisions. Simultaneous 
collisions are commutative. 



Using the simplified expression, the cluster-cluster col- 
lision is calculated in the following process. The initial 
wavefunction is simplified as 1 1 • • • It ■ ■ ' t ) with n I's 
and n f 's. The cluster-cluster collision is described by 
a series of one-to-one collisions, and (i) at the collision 
between I and t , the wavefunction splits into reflec- 
tion term (amplitude p) and transmission term (ampli- 
tude r), (ii) at the collision between two particles of the 
same spin, the wavefunction is multiplied by —1. The 
ordering of the one-to-one collisions is the same as in the 
classical dynamics, shown in Fig. |4l The first collision 
occurs at the center of the system, and the second colli- 
sions occur at the locations next to the center (n > 2). 
The third collisions arise at the center and at the two 
locations away from the center {n > 3). Finally all the 
particles have the reversed momenta, and then collision 
dynamics finishes. The external parameters of the ICM 
are only p and r, which are of the one-to-one collision. 
Since the calculations for large n are complicated, we nu- 
merically compute the ICM reflectance R]^^ (u) and the 
ICM transmittance T^^^(w). 

We plot the reflectance obtained by the ICM R]^^iu) 
in FiglDJa) and the ratio R}^^ (u) / Rf^^ (u) in FigHJb) 
for < u/uc < 1. The figure exhibits that R]^^iu) 
agrees fairly well with i?„(u) (the DMRG simulation loses 
accuracy in the limit of u ^ 0), and the ICM is valid 
in the region. We discuss the asymptotic behavior of 
R]^^{u) as m — 0. In this limit, the dominant terms 
of R^^{u) are the coefficients of single- reflection wave- 
functions (wavefunctions after only one reflection), since 
almost all particles transmit in this limit. The coefficients 
of the single-reflection wavefunctions are pr""^, and the 
number of the single- reflection wavefunctions is n^, so 
RI^^{u) ~ n^|pr"-^p/n ~ n{u/uc)^ is the asymptotic 
form. This analytical calculation is free from the interfer- 
ence effects, so the absence of quantum effects in the limit 
of M — )■ is supported by the ICM. Fig|3] also shows the 



ICM resuh TI™{u) (a) and the ratio TI™{u)/TI^^{u) 
(b) for 1 < u/uc < 2^. The figure shows that T^'^^iu) 
is consistent with T„(u) (the difference comes from the 
momentum distribution and the uncertainty in the loca- 
tion of the collision), so the ICM is correct in the region 
too. The asymptotic behavior of T^^{u) in the limit 
of u — > oo is dominated by the coefficient of the com- 
ponent with a single transmission, |J, ••• itit ■■■ t)- 

h p^"^^), so in the 



^" ^)P/'n. ~ n{u/uc) ^ is the asymptotic form. 



The coefficient is p" ^r(l + p^ -f • • • 
limit of u oo for fixed n, T}^^{u) 

Different from the small u case, this term contains the 
interference effect. In the semiclassical picture, in which 
all the horizontally aligned segments of the particle tra- 
jectories in FigHlare independently spin-down or spin- up 
at some probability determined by the previous stage, 
the interference vanishes and the transmittance becomes 

|p"-V|2(|l|2+|p2|2^...+ |^2«-2|2)/^^ ThuS 

the quantum transmittance is n times larger than the 
classical transmittance, so the ICM unveils the quantum 
effects in the large u limit. 

The ICM is consistent with the DMRG simulation, but 
we have assumed the independent ordered collision to de- 
rive the ICM. To support this assumption, we have also 
studied the following system. Initially n fermions per 
spin are trapped separately, and the parameters of the 
particles are the same as those in the DMRG simulation 
of cluster-cluster collision (FiglSj. However, the trap po- 
tentials of both spins are not harmonic, but the combina- 
tion of two harmonic potentials and the distance between 
these harmonic potentials is d. The harmonic centers of 
the spin-down (up) particles are at x = —D, —D — d 
{x = +D,+D + d) (D = 5r]). 
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FIG. 5. (Color online) Initial particle density and trap poten- 
tial at n = 2 and d — Arj. 



At i = we suddenly give a momentum p to the 
spin-down particles and a momentum —p to the spin- 
up particles {p = mujD), and simultaneously we switch 
off the trap potentials and set the contact interaction be- 
tween spin- up and spin-down fermions as uS{xi — X2). 
The two clusters start moving at the velocity — p/m 
and they collide at a; = 0, then we calculate the re- 
flectance i?^"''(M, d). Except for the difference in the 
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distribution of momentum or the decay of wavefunc- 
tion, i?5j"''(u,0) should reproduce To simulate 

the system, we use the same method under the same 
condition as in the DMRG simulation above, except 
the following conditions. We calculate the time evolu- 
tion up to t = {2D + d)/{v — 2\/mujh/m) for n = 2, 
u/uc = 0.5,1.0,2.0 and < d < 4?7 [s/ mujh is the ini- 
tial momentum variation). We take 299 sites numbered 
-149, -148, • • • , -1-148, -1-149 at regular intervals, but the 
lattice constant 6x = O.ljy is the same. The result of 
the simulation is that - i?^'"*(u,0)| is under 

1.5% for all u, so we conclude that the initial particle 
distribution has little effect on collision dynamics. Thus 
if we calculate the cluster-cluster collision, we can sim- 
plify the system by assuming that initially the fermions 
are independently localized and independently collide. In 
the simplified dynamics, the wavefunction of the system 
splits after every independent collision, and the ordering 
of the independent collisions is the same as in the classical 
dynamics. Therefore the assumption of the independent 
ordered collision is numerically supported. 

As an application of the ICM, we discuss the effects 
of the interaction sign in the cluster-cluster collision. In 
one-dimensional collision dynamics, the sign of the con- 
tact interaction u does not affect the dynamics Q , and 
this property can be explained by the ICM. The imagi- 
nary unit i is contained only in p and r in the ICM. Since 
physical quantities such as particle density are real, the 
values of these quantities are not changed if we substitute 
— i for i. After the substitution, p{u) and t{u) become 
p{—u) and r(— u). This implies that the physical quan- 
tities do not depend on the sign of the interaction. 

Another application of the ICM is the repeated cluster- 
cluster collision under the harmonic potential, in which 
the fcth collision occurs at (fc - l)T/2 < t < kT/2. The 
dynamics of the multiple cluster-cluster collisions was 
simulated 0; the motion of center of mass is a linear 
function of < in a region of strong interaction. We simu- 
late this multiple cluster-cluster collision by the ICM, 
and calculate the wavefunction after the fcth collision 
by using the final state after the (fc — l)th collision as 
the initial state. We use their system parameters 0, 
and assume the location of the localized wavefunctions 



X-., = -D-{i-{N+l)/2)r]&ndx^ = D + {i-{N+l)/2)r] 
{i > 0). Thus we obtain the results that agree very well 
with their Fig. 3(b). Therefore the ICM is also useful for 
simulating the dynamics of the multiple cluster-cluster 
collisions. 

In summary, using the time-dependent density matrix 
renormalization group method and the Fermi-Hubbard 
model, we have calculated the collision dynamics be- 
tween two fermion clusters with the contact interaction. 
We have introduced the independent collision model for 
cluster-cluster collision, and theoretically checked the va- 
lidity of our model by numerically showing that the ini- 
tial distribution of particles does not affect the collision 
dynamics. Our model has reproduced the simulation re- 
sults, and explained the large enhancement in the trans- 
mittance at strong interaction. Furthermore, we have 
demonstrated its potential applications to the interaction 
sign effect and the repeated collision dynamics. 
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